A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials. 相似文献
The authors report zircon U-Pb geochronological,whole-rock geochemical and zircon Lu-Hf isotope data for the hornblende gabbro within the Khanka Massif,with the aim of constraining its formation time and petrogenesis. The zircon U-Pb dating shows that ~(206)Pb/~(238)Pb ages of zircons from the hornblende gabbro range from 120 to 129 Ma,yielding a weighted mean age of 123 ± 2 Ma,i. e.,the Early Cretaceous. The hornblende gabbro has SiO_2 of 44. 77%--46. 58% and belongs to the tholeiitic series on FeO~t/MgO-SiO_2 diagram. It displays a right-inclined REE pattern with( La/Yb)_N ratios of 3. 44 to 4. 42. The trace element spidergram shows that they are enriched in large ion lithophile elements( LILE) such as Rb,Th,U,K and Pb,and depleted in high field strength elements( HFSE) such as Nb,Ta,Ti and P,indicating an affinity to arc igneous rocks.The ε_(Hf)( t) values of zircons vary from -2. 6 to + 3. 9 and Hf model ages( T_(DM1)) range from 622 to 883 Ma.These geochemical characteristics indicate that primary magma of the hornblende gabbro could be derived from partial melting of young mantle material accreted during the Neoproterozoic. Combined with the Early Cretaceous igneous rock assemblages in NE Asia. It is concluded that the hornblende gabbro formed in an active continental margin related to the westward subduction of the Paleo-Pacific Plate beneath the Khanka Massif. 相似文献
In many arid ecosystems, vegetation frequently occurs in high-cover patches interspersed in a matrix of low plant cover. However, theoretical explanations for shrub patch pattern dynamics along climate gradients remain unclear on a large scale. This context aimed to assess the variance of the Reaumuria soongorica patch structure along the precipitation gradient and the factors that affect patch structure formation in the middle and lower Heihe River Basin (HRB). Field investigations on vegetation patterns and heterogeneity in soil properties were conducted during 2014 and 2015. The results showed that patch height, size and plant-to-patch distance were smaller in high precipitation habitats than in low precipitation sites. Climate, soil and vegetation explained 82.5% of the variance in patch structure. Spatially, R. soongorica shifted from a clumped to a random pattern on the landscape towards the MAP gradient, and heterogeneity in the surface soil properties (the ratio of biological soil crust (BSC) to bare gravels (BG)) determined the R. soongorica population distribution pattern in the middle and lower HRB. A conceptual model, which integrated water availability and plant facilitation and competition effects, was revealed that R. soongorica changed from a flexible water use strategy in high precipitation regions to a consistent water use strategy in low precipitation areas. Our study provides a comprehensive quantification of the variance in shrub patch structure along a precipitation gradient and may improve our understanding of vegetation pattern dynamics in the Gobi Desert under future climate change.